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Abstract 

A mechanism has been suggested recently to generate the neutrino mass out of a dimension- 
seven operator. This is expected to relieve the tension between the occurrence of a tiny neutrino 
mass and the observability of other physics effects beyond it. Such a mechanism would in- 
evitably entail lepton flavor violating effects. We study in this work the radiative and purely 
leptonic transitions of the light charged leptons. In so doing we make a systematic analysis of 
the flavor structure by providing a convenient parametrization of the mass matrices in terms of 
independent physical parameters and diagonalizing them explicitly. We illustrate our numeri- 
cal results by sampling over two CP phases and one Yukawa coupling which are the essential 
parameters in addition to the heavy lepton mass. We find that with the stringent constraints 
coming from the muon decays and the muon-electron conversion in nuclei taken into account 
the decays of the tau lepton are severely suppressed in the majority of parameter space. There 
exist, however, small regions in which some tau decays can reach a level that is about 2 orders 
of magnitude below their current bounds. 
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1 Introduction 



The tiny neutrino mass and significant lepton mixing can be incorporated in the three canonical 
seesaw mechanisms 01 HI HI. From the point of view of effective field theories they correspond 
to the three possible realizations at the tree level [4] of the unique dimension-five operator that 
induces a neutrino mass Q. The tininess of the neutrino mass is generally attributed to the 
existence of very heavy new particles or very small couplings between the new particles and 
those that we already know of. In such a circumstance it is usually hard to detect other effects 
beyond the neutrino mass. 

The above tension between the occurrence of a tiny neutrino mass and the testability of other 
physical phenomena can be alleviated by postponing the appearance of higher-dimensional op- 
erators relevant to the neutrino mass. There are two basic approaches to accomplish this. One 
can compose new fields so that the operators first occur at one J6J, two [d 0, or even three 
loop order 0. Since the loop effects provide additional suppressing factors besides a product 
of multiple couplings, one may gain in the couplings between the new and known particles. In 
the second approach, one introduces several new fields that belong to certain high dimensional 
representations of the gauge group. To induce an effective mass operator one has to go through 
several steps to connect those fields to the light lepton fields which are in low dimensional rep- 
resentations. In this multistep seesaw, a tiny neutrino mass can be induced without requiring all 
new particles to be very heavy or their couplings to light particles to be all small. 

A realistic model in the second approach has been recently proposed in Ref [10]. It intro- 
duces a vectorlike fermion triplet and a scalar four-plet so that the effective operator responsible 
for a neutrino mass first appears at dimension seven. The potential signatures of the new parti- 
cles at the Tevatron and LHC have been studied with a special focus on the leptonic decays of 
the triply charged scalars. The idea of employing higher-dimensional representations has been 
further pursued in Ref |[TT|. where the neutrino mass is induced from a dimension nine opera- 
tor. For a systematic effective field theory approach to neutrino mass operators of a dimension 
higher than five, we refer to Ref lfT2l . 

Any mechanism for the generation of a neutrino mass and mixing is necessarily correlated 
with the physics of charged leptons. Before one can be sure that the physical processes involv- 
ing new heavy particles at high energy colliders are relevant to neutrino physics, it is necessary 
to examine that the parameter regions assumed in the analysis of high energy processes are re- 
spected by precision low energy tests in the charged leptons. Particularly relevant in this respect 
are lepton flavor violating (LFV) decays of charged leptons and muon-electron (/le) conversion 
in nuclei that are severely suppressed in the standard model (SM). The experimental bounds on 
LFV decays of the muon are already very stringent lfT3l[l4"1| . and the sensitivity to its radiative 
decay is expected to be upgraded by orders of magnitude in the MEG experiment within the 
next few years lfT5l . Significant progress has also been made in LFV decays of the tau lepton, 
thanks to the large data sample collected in recent years at the B factories [[T6l [TVl |T8l [T9J . Asso- 
ciated with the radiative decays of charged leptons are the precise measurements of or stringent 
bounds on their electromagnetic dipole moments ll2"0l I2T1 |22"1| . For \le conversion in nuclei, 
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the current most stringent constraints arise for titanium and gold l|23l [241 . PRISM/PRIME is 
expected to enhance their experimental sensitivity by several orders of magnitude in the fu- 
ture [25] . These bounds will provide strong constraints on the parameter space that will be 
useful in assessing the feasibility of detecting collider processes relevant to the neutrino mass 
generation. This motivates us to do a systematic investigation of the LFV transitions of the 
charged leptons in the model of high dimensional representations ifTOll .LFV decays have been 
previously studied in a similar fashion in various models of neutrino mass generation, like su- 
persymmetric models 11261 . seesaw models Il2~7ll28ll2~9ll . mirror fermions ll30ll . little Higgs OTl . 
and color-octet particles ll32l . to mention a few amongst many. Reader should consult Refs. 
Il33l l34l for a more complete list of literature on the subject. Similarly, the \le conversion in 
nuclei has also been widely considered in many scenarios of new physics beyond SM, such as 
supersymmetric models [|35l . seesaw models ll36l . littlest Higgs with T parity 11371 . Z' models 
ll38l . and so on. The formulas and elaborate discussion of \ie conversion in nuclei have been 
given in Refs. [13311391 l40l. 

In the next section we shall make a complete analysis on the flavor structure in the model 
of high dimensional representation. The mass matrices are parametrized in terms of physical 
parameters and then diagonalized approximately. The LFV decays, the contribution to dipole 
moments of charged leptons, and \ie conversion in nuclei are then calculated in Sec(3j In Sec. 
|4] we illustrate our numerical results by sampling over a few parameters that are potentially 
interesting. We discuss and conclude in the last section. 



2 Model 



To avoid the occurrence of the dimension-five neutrino mass operator at tree level, one should 
exclude the fields that carry the same quantum numbers as those in the three canonical seesaw 
models. Since the neutrinos are in the doublet representation, the easiest approach would be to 
arrange a Yukawa coupling that connects the neutrinos to a new scalar field and a new fermion 
field which differ in weak isospin by 1/2. One way to accomplish this is to introduce a scalar 
multiplet with weak isospin 3/2 and a fermion multiplet with weak isospin 1. This avoids the 
type 1 and type 2 seesaws automatically, while the type 3 is avoided by assigning a different 
hypercharge to the fermion multiplet. The is indeed the basic idea behind the model building in 
Ref. ifTOll. The new fields are denoted as 



4> 



/ *3 \ 
4>2 



(3/2,3/2), E 




:u), 



(i) 



where the numbers in the parentheses stand for the weak isospin / and hypercharge Y /2, respec- 
tively, and the subscripts to the fields indicate the electric charges in units of \e\. The fermion 
fields E are assumed to be vectorlike to avoid chiral anomaly. The relevant SM fields are the 
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Higgs doublet and the lepton fields (with the subscripts L, R denoting chirality), 

H= {h°) t 1 / 2 ' 1 / 2 )' Fl= {7l) (1/2 '" 1/2) ' fR (o '- i) - (2) 
2.1 Yukawa couplings and mass matrices 

The neutrality in the hypercharge allows the following terms and their Hermitian conjugates: 
Flf R H, F L f R <i>, F L LH*, F L E*4>, £*E. It is possible to assign the lepton number, L(f R ) = 
L(Fi) = 1, L(E) = — 1, L(<t>) = —2. Then L is violated when 4> develops a vacuum expectation 
value (VEV). Now we write the terms in a form that respects SU (2)l and Lorentz symmetries. 
The first and last terms are trivial, Fl/rH, EE, while the second one is forbidden. For the third 
term, H* should be replaced by H = i<72H* to preserve its identity as a doublet. To form a 
Lorentz scalar without complex conjugation out of two spinor fields, we can use the charge- 
conjugated fields. The required invariant form for the third term is, in terms of Clebsch-Gordan 
coefficients, 

(^ E )o = ^(-^ E2+ ? //0 * E °)-7f(^- ; 5^) Ei ' (3) 

where the subscript on the left indicates its weak isospin of the product. We have used the 
notation for charge conjugation that = ( Vl) c and y/- = ^y°i//* with ^ = ?'y°y 2 . 

For the fourth term in the list we note that, since the vector representation of SU (2) is strictly 
real, (E^E^Ej) is a vector when (E2,Ei,Eo) is. The invariant form is thus, 

(eX^ = I(E^n L 4> -Ei/L4> 3 )+^(^L^2-WL*l) 

+-|=(ET/l^2-^ l 4> 1 ). (4) 

Including the generation index in SM, the mass terms and Yukawa couplings are summarized 
as follows: 



- ^Yuk+mass = m E EE + [v, ,/•}.,/ A>/ // + Xj (*gi5x) q + Zj (zF Lj <S>) o + h.c. 



(5) 



where y is a 3 x 3 complex matrix, x, z are each a three-component complex column vector, and 
m£ is real positive by definition. 

When the electric neutral components of H and 4> develop a VEV, 

the Yukawa terms will contribute to the masses of the neutral and singly charged fermions. For 
simplicity we shall assume in this work that both VEV's are real positive. The mass terms are 
written as 

=-^mass = =-^mass,2 ~t~ =^mass,l -^mass.O, (7) 
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with the number in the subscript denoting the electric charge. The field E 2 has only a bare mass: 



-_Sf 



mass,2 



m s E 2 E 2 , 



(8) 



while the fields of other charges also derive masses from Yukawa couplings so that mixing 
between the light and heavy particles can appear. For the singly charged fields, we have 



where the four-component column fields and the 4x4 mass matrix are 



(9) 



IR 



fR 
If 



1L 



h 
if 

^IR 



Mi 



V2 



V2 y 





2^3' 



(10) 



Since the neutral fermions are generically Majorana, their mass terms are apparently more com- 
plicated. With the help of charge-conjugated fields, they are 



~^mass,0 = 

where the fields and mass matrix are 

VOR 



(11) 











11) 


l,*«=l 




,M = 











^z T m z 



V 2^2 



(12) 



These mass matrices will be analyzed and diagonalized in the later subsection. 



2.2 Scalar potential 

Before we diagonalize the mass matrices of leptons we discuss the scalar potential for com- 
pleteness. Remember that the fields <I> and H have the quantum numbers I = Y /2 = 3/2, 1 /2, 
respectively. The possible quadratic terms are, H^H and 4>^4>, while there can be no trilinear 
terms. For the quartic terms there are two sets of them: either <3> and 4^, and H and H* come 
in pairs, or one <t>^ is accompanied by three H. For the first set, the following ones are obvious: 

(la) (H f H) 2 , (4> f 4>) 2 , (H^H) (13) 

The pure H term is unique. This can also be understood as follows: with two identical H's of 
7=1/2 one can only construct an I = 1 form that is symmetric in the two H's. (The 1 = 
form vanishes identically.) This is also the case with two H's. From the two 1=1 forms one 
can construct a unique 1 = term (that is symmetric, though not necessarily, in the two forms). 
Indeed, it is easy to check that x a HH^ x a U = (H^H) 2 where x a are the Pauli matrices. But 
this is not the case with the <I> field. With two identical 3>'s of I = 3/2 we can construct two 
forms that are symmetric in them, one of I = 3 and the other of I = 1 . This implies that there 
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are two independent, invariant, pure <I> terms. Similarly with the half-// and half-4> terms. The 
additional terms in the first set are thus 

{lb) (<P% a /2 <P)(<t>% a /2 <P), (//V//)(4>tr 3 fl /2 4>), (14) 

where Ty 2 stand for the generator matrices for / = 3/2. 

For the second set of quartic terms involving three //'s and one <£, we first combine three 
H's into a form of / = | which must be symmetric due to Bose symmetry. Then, out of this 
form and one <I> field, we form an / = term: 



(*HHH) = ^o(H°n 3 + ^iH-(H *) 2 + ^ 2 (H-) 2 H°* + ^ 3 (H-) 3 . (15) 



The normalization in the above term looks a bit unusual due to the appearance of identical fields. 
The most general scalar potential is thus 

+A// t //4> t $ + ^A / (// t T a //)($ t r 3 fl /2 $) + [k(<PHHH) Q +h.c] , (16) 
where all couplings except K are real. We note in passing that the X 1 ^ term was missing in Ref. 

ma. 

The VEVs of the scalar fields are determined by requiring the vanishing of the first deriva- 
tives and positive-definiteness of the matrix of the second derivatives of the potential. We can 
always choose, by a global U(l)y transformation, one of the VEVs, say, V2, to be real positive. 
Then one can see that KV4 must be real, and the vanishing conditions become 

13 3 
+ X H y 2 2 + 2^ M 2 + g^|v 4 | 2 + -v 2 kv 4 = 0, 

-^M*|v 4 | 2 + v 4 | 4 + |a4|v 4 | 4 + ^Av|| v 4 | 2 + ^^^^412 + ^fcv 4 vi = 0. (17) 
z z o 4 lo o 

Since V4 7^ breaks the custodial symmetry, it is natural to assume | v 4 1 V2. Assuming further 
that the quartic couplings are perturbative, we have to good precision that 

v 2 *M ^« 2 J^i (18) 

VV v 2 U H iil-(4X+3X>)n 2 H 

Since the K term breaks lepton number, it would be easy to attribute the tininess of V4 to that of 
K. With so many free parameters at hand it is no problem to guarantee that the above is the true 
vacuum. If k vanishes, it would be necessary to fine-tune the parameters to get a tiny V4, which 
we shall not pursue further. 

2.3 Diagonalization of lepton mass matrices 

We continue to assume for the sake of simplicity that both VEVs are real positive. To diag- 
onalize the mass matrices, we first parameterize them without losing generality in terms of 
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independent physical parameters by following the procedure advocated in Ref. BTl . We sketch 
below how this is done. 

With three generations in SM there is always one massless neutral mode, which can be Vi [in 
normal hierachy (NH)] or V3 [in inverted hierarchy (IH)] according to the oscillation data. (With 
n g > 3 generations there are n g — 2 massless modes while there is none with less generations.) 
We describe the case of NH in some detail and will record the result for IH later. By applying 
a judicious unitary transformation to nu, we can convert the column vectors x and z into the 
standard form: 

X = (0,0,x) T , Z=(0,z,c z ) T , (19) 

where x, z are real positive and c z is generally complex. This fixes the phases of the two fields 
(named again as nu_ and n.13) orthogonal to the massless mode (nu) but leaves the latter's phase 
free. To keep the partnership under SU (2)/, between the neutral n± and charged fi fields gen- 
eration by generation, the same transformation should be applied to fij as well. This modifies 
the entries in y but does not alter its generality. We assume that this has been done already. 

To reduce the y matrix to its minimal form, we proceed as follows. By a unitary transforma- 
tion of fRj, we can cast y into the form: 



(20) 



where V2 and V3 are real positive and C2 is complex. This fixes the phases of ^3,^2 but leaves free 
that of fin. By rephasing further fn — > e'Pfu,, which is augmented with nu — > e l ^nn to pre- 
serve the SU (2)i partnership, the first row of y becomes effectively, e~ l $ (jie mx , r2e l(X2 ,r?,e ia3 ). 
Now we choose /3 = «2 (or equally well, /3 = (Z3), which fixes the phase of fn and thus that 
of the massless nn as well, and then choose <Xi = /3, which fixes the phase of fm. This leaves 
with us the final version of the y matrix: 





f rie iC(i 


r 2 e iai 


r 3 e ia i 







yi 






V 





3'3 



y\ 


)'4 


Cl 





yi 


Cl 








3'3 



(21) 



where ^1,2,3,4 are real positive and ci,2 are generally complex. 

We have used up the degrees of freedom in defining fermion fields to reduce the num- 
ber of parameters in the Yukawa sector to its minimum, i.e., in terms of independent physical 
parameters, without sacrificing generality. There are seven real positive parameters (m£ and 
x i 71,2,3,4) z) and three complex ones (ci,2, z )- They will be traded for nine masses (of one 
doubly charged, four singly charged, and four neutral fermions), three CP phases, and a single 
independent mixing angle. It looks challenging at first sight for the model to accommodate the 
two large mixing angles measured in oscillation. But as we shall see later, all heavy fermions 
are nearly degenerate, which effectively saves parameters at our disposal. The mass matrices 
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for the singly charged and neutral fermions are summarized as 

/ 3 



Mi = mz 



£2 v 
0lx3 



£2 
2%/3 
1 



X 



, Mq = mz 



7S xT 1 



(22) 



where £2 = V2/mx, £4 = v^/mz with £4 being necessarily tiny. The submatrices X, Y 1 Z for NH 
have been given in Eqs. (IT9T ) and (12TT) . while for IH they are 






(23) 



Now we diagonalize the mass matrices Mq and M\ perturbatively exploiting the hierarchies 
in parameters. We describe how this is done for the NH case and will indicate how to obtain 
the analogous results for the IH case. To save writing, we put a prime to all parameters in the 
standardized Yukawa matrices shown in Eqs. (fT9b . (1271) . and d23l . and reserve the unprimed 
parameters to those that have been multiplied by a factor of £2,4; namely 



e 2 / / v 



^=(z',4), {y h a) = ^iy'jA)- 



(24) 



The symmetric complex matrix Mq is diagonalized to real nonnegative eigenvalues by a unitary 
matrix U : 



U MqU = diag(mi,m2,m3,m4,m5^ 



(25) 



where m\ = for NH while = for IH. U is composed of the five column vectors u(j) 
corresponding to the eigenvalues mj, U = («(1), • • • ,m(5)). Since £4 is tiny, we can solve the 
problem perturbatively in the z, c z parameters (collectively denoted as £) while keeping exact 
dependence on x. After some work, we obtain the eigenvalues for NH: 



^ 0, 

: so 

— = Sq 



m\ 

m 2 

m 3 

m,4 

mz 
m 5 

mz 



Z 2 + ci\c 



-co c« 



z 2 + cl\c z \ 2 + c \c z 



+ 0(e 3 ), 

+ 0(£ 3 ), 





1 


\C Z \CQSQ + 


2 C0 




1 


\c z \cqSq + 


2 C0 



cIz 2 + {\-34s Z q)\c z \ 1 +6>(£^) : 



,2„2 



,2_2 



2„2\ 



(26) 
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corresponding to the light eigenvectors 
u{\) - 



;i,o,o,o,o) r , 
/o 



u{2) 



s m (-c e 



-ia- 



-co\Jz 2 +cl\cz\ 

-ia. 



s 0\ c z 



u(3) 



and the heavy ones 



-ia- 



\ s m s Q e~ 

/o 

c m cae~ 

C m CQ ( - CqsJz 2 +cl\c z \ 2 +sl\c z 

V c m (-s e- ia z) 



/0 



+ 0(e 2 ), 



+ 0(e 2 ), 



»(4) 



\ 



c 2 z 



s -2C§(2cg-jg)|c z 



1 + 5Cq5o|c 



c - 5Coso(2so-Co)k z | 



+ 0(e 2 ), 



/0 



»(5) 



\ 



c 2 z 



5 +2 C o( 2c O- S o)l^ 



1 - \cIsq\c z \ 



\ 



co + \cqsq{2sI-cI)\c z \ 



+ 0(e 2 ). 



J 



(27) 



(28) 



Here we have parametrized c z = \c z \e' a: and e'^ az+2 ^ = ±1. And the triangular functions are 
defined as follows: 



1 



so 



Vl+x 2 



CO 



: 5 ^» 



m2 



m 3 



One can tidy up the 0(e) terms in u(2) and w(3) in terms of mass ratios using 



(29) 



1 1 



co y z 2 + c 2 , I c z 1 2 + ( -) s 1 c z I J = — — [m 3(2 ) + (c - s ) m 2(3) ] 



(30) 



but no similar simplification occurs for the terms in w(4) and w(5). 

We make a few comments on the above result. The light neutrinos gain a mass of order 
m ~ x'(z' or |c^|)v2V4/m£. For a scalar potential with /ij, 3> /i^, corresponding to 4> particles 
much heavier than the SM Higgs, we have from the previous subsection that V4 ~ \ k\v\/ ji^. 
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This yields a neutrino mass that is triply suppressed by heavy scales as designed in Ref. IfTOl . 
Second, the two heavy neutrinos are almost degenerate with a splitting as small as the one in 
light neutrinos: 

»?5 — tri4 pa «?3 — rri2 ~ 2co5ol c zl m z- (31) 

The significance of this at high energy colliders will be explored in the future work. Third, 
the result for the IH case whose Yukawa coupling matrices are parametrized as in Eq. (l23~l) is 
obtained by first reshuffling the labels of the light solutions, {m\ 1 m'i 1 m-$) — y (m^, mi, m 2 ) and 
(w(1),m(2),m(3)) — > (w(3),w(l),w(2)), in Eqs. (12"6l) and (|27T) . and then interchanging the first 
and third rows in all u(j). 

Now we diagonalize the mass matrix Mi for the singly charged leptons by bi-unitary trans- 
formations, U^ l M\Ur = diag(m e ,m jU ,m T ,m^), where % is the heavy lepton of charge — \e\. Since 
the matrix y is not related to new physics, its entries should be naturally small. At this stage 
we have no idea on how large the parameter x could be, and thus leave it free. We therefore 
diagonalize M\ in two scenarios according to whether the x parameter (scenario A) or the y 
parameters (B) are treated perturbatively. As will be clear later, the x parameter is severely 
constrained by LFV transitions of the muon so that both scenarios serve as almost equally good 
approximations. We shall describe the diagonalization for the NH case and indicate at the end 
how to obtain the results for the IH case. 

In scenario A we do perturbation in x. We first solve the zeroth order eigenvalue problems 

(u L (e) u L (ii) u L (x))^l(u L (e) u L (fi) u L (z)) = diag(A„ A M , 

(u R (e) u R {v) u R {x))K^l{u R {e) u R {^i) «*(t)) = diag(A e , A M , A T ), (32) 

for the two 3x3 Hermitian matrices which share the same set of real positive eigenvalues: 

„2 i „2 , \„ |2 



yi+yi+\ci\ yiy^ J rc\c 1 y 3 ci 

,2 _ | ^ ^ J ± u \2 

''L 



c.c. J2 + F2I 

C.C. c.c. 



/ y\ y\yA y\c x 
Jtl = c.c. yl+yj yic 2 + yACi | . (33) 

\ C.C. C.C. |ci| 2 + |c 2 | 2 +J3 

We denote the light charged leptons by the beginning Greek letters a and /3, and introduce the 
auxiliary variable and vectors: 

8a = (1 - ^ L = -Y,5 p u I ?(P)u L (P), Z R = -^dpu** (j8K(j8), (34) 

P j8 
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where the subscript in Ua R {fi) denotes their a-th entry. The eigenvalues are found to be 



4 = Xa-^-\u L r (a)\ 2 8 a k a + 0(x 3 ) 
x 2 

= K- -yj\u R (a)\ 2 8 a + 0(x 3 ), a = e,\x,x (35) 

yw 2 

— | = l + -^\u L T (a)\ 2 8 a + 0(x 3 ) 

= l + ^ + ^L*>)| 2 5« + 0(A (36) 



2 2 



a 



where the two expressions from diagonalization are equivalent using y 2 \u R (a) \ 2 = X a \u^(a) | 2 , 
and the unitary matrices exact to 0(x) axe 

( u L (e) u L (ii) u L (x) -j=% L \ 

Ul = \ tM^) f 2 ^ L M fM^ 1 J^ 1 ' 1 ' 1 '^' (37) 

( uR{e) uR ^ ] uR{%) fin ^ ™ 

% " ^«(e) ^5 Ma f(M) ^«(r) 1 Jd iag (l,l,l,^). (38) 

In the above, the 0(x) phases pm are not fixed completely, but PrP* l is fixed by requiring that 
the eigenvalue, m x , in U^ l M\Ur be indeed real positive. The same arbitrariness also occurs in 
u L R , but the number of physical phases in Ul,r is restricted to two by the phases appearing 
in M\. Considering A« 1 in practice, we have to good precision, 8 a = 1 + 0(X a ), so that 
%y ,R = — 5 T y + O(Ace) from unitarity of u L and u R , which simplifies things a bit. 

In scenario B we instead treat the parameters (y ; ,c ; ) (denoted collectively as 7]) perturba- 
tively. By solving first the eigenvalue problems, 

(u L (e) u l (li) u L {x))^l{u L {e) u L {\i) u L (x)) = diag(A e , A M , At), 

(u R (e) u R {pi) u R {x))^i{u R {e) u r (^l) u r (x)) = diag(A e , A M , A T ), (39) 

where 



•^X — l>'3^c + >'3 5 -^j? — -^7? |y3— ^c + y3 j (40) 



we find the eigenvalues 



.2 



2 



A a + 0(77 4 ), a = e,jU,T, (41) 



m ^ i , 1 2 , ^/„2> 



l + ^ + Ofa 2 ), (42) 
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and the diagonalization unitary matrices: 
/ u L e {e) u L M 



U, 



( 



U R 



u R e {e) 
u*{e) 
•*f e ) 



c+4(t) 

M f(jU) 
K?(/i) 







(43) 






-y3J+c+ 
l 



(44) 



V y3J+c+«f (e) y 3 s + c + u*(Ll) y 3 s + c + u^(T) 
The triangular functions in this scenario are 

s + = , c + = , (45) 

When all parameters x and (yi,ci) are small in magnitude and treated on the same footing, 
both scenarios yield an identical result to the leading order. In both scenarios, the following 
mass relations among the heavy E particles hold: 



m zV 1 +x 2 . 



(46) 



m z, m X — m sy 1 +^ 2 /2, m V4 ~ mv 5 

In our later numerical analysis we shall work with scenario B. The explicit results displayed 
above are for the NH case. For the IH case whose matrices are parametrized as in Eq. (1231) and 
without changing the increasing order of the mass eigenvalues, Ulr are obtained from those 
for NH by interchanging the (1,3) rows which are computed with the parameter interchanges 
yi <H- V3 and c\ -H- C2 made. 



2.4 Couplings of leptons 

The gauge interactions of the leptons are 



(47) 



The currents are written first in terms of weak eigenstates and then grouped into M^i^il and 
cw4 = WLzlf^0L+WLzif^iL + ^ZRY^iR-(l-2si)T^r^, 



where i stands for all four leptons of charge 



1. The coupling matrices are 





3 



z\ 



(48) 



(49) 
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with the usual notation sw = sin 6w and cw = cos Oyy- In terms of mass eigenstates, ^ol = U Vl, 
^or = U*Vr = U*v£, = Ul£l, and = Ur£r, the weak currents finally become 

Cw j£ = w^ L VvL+4^ L 1 r M 4+4^ 1 r M ^-(i-24)^r M s2, m 

where 

w L ° = u^w° L u L , w£ = u T w° R u R , wl = uy L , w£ = uy R , 

2$ = Uh° L U, 2ft = ulzlu L , 2?r X = uIzrUr = s 2 w U. (51) 

The upper-left 3x3 submatrix in sflW^J (y/2 from our normalization convention for cur- 
rents) corresponds to the Pontecorvo-Maki-Nakagawa-Sakata (PMNS) matrix in the seesaw 
limit, and will be called the effective PMNS matrix or Vpmns- To help understand the procedure 
to be adopted in our later numerical analysis, we take a closer look at it. As both scenarios 
A and B are similar for a not large x parameter, we illustrate our discussion in the latter. The 
submatrix is 

Vpmns = C/J L diag(l,l,c + /c )C/ n . (52) 

Here UfL= {u L (e) « L (/i) u L {x)) is the unitary matrix that diagonalizes while U n is the 
unitary matrix that diagonalizes the light neutrino mass in the seesaw limit upon choosing 
e~ l P+ = e'" / 2 and e~'P~ = ie ia ^ 2 . In other words, we have Vpmns = uj L U n in the seesaw 
limit, and c+/cq ps 1 + jc 2 /4 measures the small departure of Vpmns fro m unitarity when the 
mixing with the heavy particles is taken into account. As the unitarity has been checked at a 
precision not better than a percentage, it is safe if x assumes a value not larger than, say, 0.1. 
This will be fully respected in our numerical analysis. This result applies to scenario A as well 
where x is small by definition. 

We briefly highlight some other features in gauge couplings relevant to LFV transitions of 
charged leptons. The light charged leptons have 0(z,c z ) suppressed couplings to heavy neu- 
trinos in #2°- The massless neutrino decouples from W R , while all light charged leptons are 
suppressed in it by a factor of xy3. A similar situation occurs in the charged current involving 
singly and doubly charged leptons: the left-handed part of the light charged leptons is sup- 
pressed by x and the right-handed by xyj,. Finally, the neutral current of the singly charged 
leptons is dominantly flavor diagonal, with an 0(x) mixing between the light and heavy parti- 
cles and an 0(x 2 ) mixing amongst the light leptons. 

As most new gauge couplings between light and heavy particles are controlled by the x 
parameter, we shall consider its largest allowed value in numerical analysis. The spectrum of 
the light neutrinos then implies that the z parameters must be extremely small. We can therefore 
focus in the Yukawa sector on the x terms. Ignoring the tiny mixture between the doublet and 
four-plet scalars that is proportional to V4/V2, the terms relevant to our study are 

-J^Vuk ^ -^ + hx. = -^([/;) 3fl (^W^ + k, (53) 
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where the summation is over all singly charged leptons and h is the physical scalar in H° = 
(h + iG°)/V2. 



3 Lepton flavor violating transitions 

3.1 Radiative decays and electromagnetic dipole moments 

A direct consequence of the neutrino mass and mixing mechanism in the last section is the 
lepton flavor violating transitions of the light charged leptons. We start with the radiative decay 
—> J, whose amplitude has the dipole form 



—— — y^up (p — q)(hjJ > L J r hRP^io^yUa 



(54) 



where p, p — q are the momenta of the initial and final leptons, q, e{q) are the photon's mo- 
mentum and polarization vector, and Plr = (1 =F 75) /2. All information on dynamics is stored 
in the form factors hiR. The decay width is 



ma0lGF -(\h L \ 2 + \h R \ 2 ) 



2 9 



7T 



(55) 



where we have ignored the mass of £p in phase space. 

For the model considered here, the form factors h^R are contributed by the Feynman dia- 
grams shown in Fig. [T] Figures (la)-(lc) involve the charged currents between the singly and 
doubly charged leptons, and between the neutral and singly charged leptons, respectively, while 
Figs. (Id) and (le) originate from the flavor-changing neutral currents (FCNC) and physical 
Higgs exchange. For the gauge-boson mediated graphs we compute in the unitary gauge. This 
is simplest but caution must be exercised to cope with a technical point concerning constant 
terms, see the last paper in Ref. ll30ll . We have carefully done the Dirac algebra in d dimensions 
before the limit d — > 4 is taken, to avoid missing certain finite terms, and find consistent results 
with that work. 

Ignoring again the terms suppressed by mp and keeping terms up to the linear order in m a , 
the form factors from each graph in Fig. Q]are 



fiL(a) 
h L {b) 
h L {c) 
h L (d) 

h L (e) 



-2JK 



K ja m a ^(rj)+W L ja m^( rj 



2W t 



V2^T2 U ^ UL ' 35 



(56) 
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where summation over the virtual lepton flavors is implied, and 

h R {a,b,c,d) = h L (a,b,c,d)\ L ^ R , h R {e) = h L (e) \l^r,3^4- (57) 

We have denoted the ratios of the masses appearing in the loops as r% = tn^/m^, rj = mj/m^r, 
s§ = m^/m|, t§ = mi/rn^, where j and 8 enumerate all virtual neutral and singly charged 
leptons, respectively. Some products of the coupling matrices in the above can be simplified 
using their explicit forms. For instance, the first term in hi(d) drops out since 3f R is diagonal; 
the matrices in the charged current involving the doubly charged lepton are ^ R ( L \ a = — Ur(j\ 4«- 
And the loop functions are 

&(r) = , 1 -7 [l0-43r + 78r 2 -49r 3 + 4r 4 + 18r 3 lnr1, 
6(1— r) 4 L J 

#(r) = / 1 x J -4 + 15r-12r 2 + r 3 + 6r 2 lnr1, 
(1 — ry L J 

jt?( r ) = 1 [-8 + 38r-39r 2 +14r 3 -5r 4 +18r 2 lnr], 
J{f) = ^A^[4-3r-r 3 + 6rlnr], 

= 12(l-r) 4 ^ 2 + 3r ~ 6r2 + r3 + 6rlnr ]' 
^ W = ^Y^7)3 [ - 3 + 4r - r 2 - 2 In r] . (58) 

Related to the above radiative transition amplitudes are the anomalous magnetic moments 
and electric dipole moments of the singly charged light leptons. They are worked out to the 
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linear order in the mass m a of the considered lepton. The anomalous magnetic moment, defined 

as a = (g — 2)/2, is 



a(£a] 



2\f2GFin a 
(4tt) 2 



h(a) + h(b) +h(c) + h(d) + h(e) 



(59) 



where 



h(a) 

h(b) 

h{c) 
h(d) 

h(e) 



J2 



I U L>3S U Rj4a 1 2 + I U RAS U L)3a 1 2 ) m a X (t s ) 



V2G F m 2 h 12 
+Re(Ul 4a U^ s U* s U L , 3a )m s 3?(t s ) 



(60) 



Again some of the above can be simplified using explicit forms of the coupling matrices given 
in the last section. 

The electric dipole moment of the fermion \y is defined as =^ e dm = — id/2\j/a^ v y5 i//F^ v , 
with F^ v being the electromagnetic field tensor. It is evaluated to be 



where 



d(£ a ) 



h(a) 
h{b) 
h{c) 
h(d) 

h(e) 



2V2eG f 
[An) 1 



h(a) +h(b)+h(c)+h(d)+h(e) 



(61) 



-2Im(^ 2 a #^)m^(r z ), 
2\m{W R %Wl ja )m0{r j ), 
2lm{Wi a W%)m^J?{r^), 



V2G F m 2 h 12 



(62) 



Since 3f R is diagonal and 3?l is Hermitian, there is actually no contribution from the FCNC 
graph, h(d) = 0. 

3.2 Purely leptonic transitions 

Now we consider the purely leptonic transitions of the light charged leptons. These include the 
experimentally well-bounded decays £§ — > £ a £p£y and the muon-electron (/ie) conversion in 
nuclei. The leading contributions in the model considered here arise from FCNC couplings of 
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Figure 2: Feynman diagrams for £ a — > £p£y£p- 



the Z boson. The Higgs exchange is suppressed by additional factors of x and a heavier Higgs 
mass, while we have verified that the photonic contribution is always subdominant. We do not 
consider LFV decays of the Z boson since they are experimentally much less constrained. 

There are three types of decays corresponding to a = /3 = 7, a = 7^/3, and a — /3 ^ y. 
Since the flavor-changing couplings carry a factor of x 2 , only the transitions of the first two types 
are important while the third one is severely suppressed. We thus concentrate on the decays, 
£a ^p^ylp with 7 = j8 or 7 ^ j8, whose leading terms come from diagrams shown in Fig. [21 
Note that there is a relative minus sign between the two graphs and that for 7 = /3 one should 
attach a factor of 1/2 in the total decay rate. Ignoring the final-state masses, the rate is given by 



4r 



1 



I ayl <jW \2 ■ 1 av>\ av>\ 1 
\° £ L,ya- c L,pp\ ~+~ \^L,Ba^L,yB\ 



1 + 5p y 1 

T^e\-£L,Ya° £ L,PI3° £ L,Ba° £ L,yB) 

+ 1 3?L,ya ^R,PP \ 2 + \^L,l5a ^R,yP \ 2 + (L ^ R) 



(63) 



with To = G 2 7 m 5 a /( 192tt 3 ) being the decay rate for the dominant decay mode, £ a — > V a £pVp- 
Since is diagonal and a^j3, a ^ 7, the term (L R) actually drops out. 

A competitive process is the coherent fie conversion in nuclei, /l~N — > e~N. It involves 
various atomic and nuclear effects in addition to the short-distance physics of lepton flavor 
violation. A comprehensive study has been given in Ref. ||39l based on the method developed 
in Ref. [|40l . which improved over earlier efforts on various corrections Il42ll43l to the original 
calculations Il44ll45l . These corrections turn out to be particularly important for heavy nuclei. 

The effective Lagrangian relevant to the coherent conversion via the leptonic FCNC cou- 
plings of the Z boson can be written as lT3~9l 



Gf r 

7^ £ [(gLV( q )efPLfi+gRV( q )efPR^)w 

V ^ q=u,d,s 



+ h.c. 



(64) 



where g LV{ll) = (2 - 16^/3) J^ eM , g L v(d, s ) = (-2 + 8^/3) J^ eM , and g RV{q) = for the con- 
sidered model. Then, the \ie conversion rate is 



T(fi~N e~N) = 2G 2 F 



1 \/ W _i_ ^ 1/ (") 1 2 _i_ I 7^ V ( pS > _1_ 7^ 1/ W 1 2 
\SlV v N ~ t 6LV v N I ~i~\SRV v N ~T~5RV V N I 



(65) 
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where g§ = 2g LV{u) +g LV{d) , g§ = g LV ( u )+2g LV{d) , and similarly for g$ mdg$. vjf ] and 
(n) 

' are overlap integrals of the /i, e with the protons and neutrons in the nucleus N, which have 
been numerically evaluated and cataloged in ||39ll . The above rate is usually normalized to the 
corresponding ordinary muon capture rate, Ct) ca pt, to yield a branching ratio, Br(/i~/V — > e~N), 
for the fie conversion on a particular nucleus N. Since the purely leptonic decays and the \le 
conversion in nuclei originate from the same FCNC couplings, the ratio of their branching ratios 
has the simple form, 



Br^-N^e-N) _ G 2 [(2 - Ss 2 w )vj p) - 2V^] 2 
Br(/x eee) © capt (l - As 2 w + 6s^) 



(66) 



upon ignoring minor corrections in the diagonal element 3£' L 1 . Namely, the relative importance 
of the two transitions rests on that of their experimental bounds. 



4 Numerical analysis 

As we shall see later, the Yukawa coupling x' , or x = x'\2/ (v^wz), that couples the light and 
heavy fermions via the Higgs doublet is a central parameter that controls the overall scale of the 
LFV transition rates. We mentioned earlier that the parameter measures the unitarity violation in 
the effective PMNS matrix. Since the heavy fermions have a squared mass splitting proportional 
to x 2 [see Eq.(l46l)l. it could also be sensitive to the violation of custodial symmetry measured 
by the parameter Ap = m 2 ^ / (c^m|) — 1. We have calculated the one-loop contribution due to 
the heavy fermions 

Ay/2G F ml\9 4 c /4 lTeV 2 
APe = — — -j^— r«1.7xl0 x ^ — . (67) 



{Ail) 2 48 



This is balanced by that of a nonvanishing VEV of the nondoublet scalar field, V4 7^ 0, that 
occurs already at the tree level, Ap<$ w — 6v|/v 2 (see also Ref. ifTOlO . We noted in the above 
that half of the power in x comes from the mass splitting of the heavy fermions. The other half 
originates from the mixing effect in the two vertices, which is essential for a contribution to Ap 
since a vectorlike multiplet cannot contribute even if it is not degenerate. Since the p parameter 
is measured at a precision not better than 10~ 4 , we are on the safe side if x' is not larger than 
0.7 even for a doubly charged fermion as light as 200 GeV. This is a much weaker constraint 
than we shall get below from LFV transitions. 

Before we show our numerical results, we outline how the free parameters are manipulated 
based upon the formulas in Sec. [3l Since the unitarity of the PMNS matrix has been verified to 
certain level, we can use it as a guide in browsing the parameter space. In this way we can cover 
the majority of the parameter space that is consistent with an almost unitary effective PMNS 
matrix, Vpmns- The PMNS matrix generally contains three mixing angles, one Dirac phase, and 
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two Majorana phases in the standard form: 

/ C12C13 5 i2 Ci3 StfU* s \ 

^PMNS = I -C23S12-S13S23C12U8 C23C12 ~ *13*23*12"5 *23Q3 diag(«i , U 2 , U 3 ) , (68) 
\ S23S12-S13C23C12W5 -523^12 -*13C23*12"5 C23C13 / 

where c,j = cosfyy, s,j = sin0,^, m; = exp(za ; /2), and u§ = exp(iS). We use the measured 
values for those angles (while choosing some values for the phases which are not yet mea- 
sured) and for the light neutrino masses mi 2,3 in either NH or IH. Then, the heavy neu- 
trino masses and the parameters z, \c z \ are uniquely fixed once the parameters (x,m£) 
are assigned a value. The diagonalization matrix U for the neutrinos is also fixed up to the 
phase of c z = \c z \e lClz . In particular, the matrix U n that is formed from the first three rows 
of the vectors in Eq. (1271) (for the NH case and similarly for IH) is fixed up to a diagonal 
phase matrix, U z = diag(l, e' a ^ 2 , e^ m ^ 2 \ multiplied from the left. Applying the definition 
Vpmns = UfjUn and the first equation in Eq. ( 1391 (for scenario B and similarly for scenario 

A), we have Uj^£lU z = (f/^L^VpMNsdiag^, A^, A T )Vp MNS (f/^f/ z ) where the right-hand side 
is completely known with the additional input of the light charged lepton masses. We observe 
that this procedure then determines uniquely all of the parameters 3^1,2,3,4, ci,2 as well as a z on 
the left-hand side. The diagonalization matrix Ui is thus fixed. Once the Y matrix is known, we 
can follow the formalism in Sec. [3] to find the other diagonalization matrix Ur. In a final step 
we go back to check that the effective PMNS matrix Vpmns obtained from the matrices U and 
Ui determined above does not violate unitarity beyond the allowed level. 

To get some feel about the branching ratios for the LFV transitions, we start with the results 
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Figure 4: Branching ratios as a function of x for IH in scenario B and tribimaximal mixing. 



for the simplified case of tribimaximal neutrino mixing with all phases set to zero. We assume 
mx = 200 GeV in the following discussions so that the heavy fermions are within the reach of 
LHC. In the upper panel of Fig. [3] we show the branching ratios for the muon decays and jle 
conversion in nuclei ^ 7 Au and ^fTi as a function of the x parameter for the NH case in scenario 
B, together with the current experimental bounds on them (horizontal lines). The lower panel 
depicts the branching ratios for the two x decays, X — > /i/i/i, /lee, for the same range of x, 
while other decays are severely suppressed. We see that the \le conversion on the heavy gold 
nucleus sets the most stringent constraint on the x parameter though it also inherits the largest 
uncertainty from nuclear physics. The bound from the conversion on titanium is comparable 
to that from the purely leptonic decay, ji — > 3e. For such a small x the deviation from the SM 
values of the anomalous magnetic moments in Eq. (l59l) and the contribution to the electric 
dipole moments in Eq. (I6TT) are too small to be relevant. In Fig. |4]we depict the corresponding 
results for the IH case in scenario B. Note that in this case the dominant LFV x decays change to 
eee and e\x\x. Generally speaking, when assuming tribimaximal neutrino mixing, the stringent 
bounds on the muon imply that the tau decays are not likely to be observable in the near future 
for the majority of the parameter space. 

We thus ask if there is a region in the parameter space where the muon decays and \ie 
conversion in nuclei are significantly suppressed while the tau decays are not much below the 
current bounds. To help identify the interested region it is useful to work with the leading terms 
in the limits of infinite virtual heavy masses and vanishing virtual light masses. An inspection 
of Eqs. (l56l) and (l63l) augmented with the coupling matrices displayed in Sec. [2]tells us that in 
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Figure 6: Branching ratios as a function of 0:2 for IH in scenario B and best-fit mixing angles. 
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Figure 7: Branching ratios sampled over (x, OC2) (upper panels) and (x, 8) (lower panels) for NH 
in scenario B and best-fit mixing angles. 

the NH case the radiative decay £ a — > in y is dominated by the terms with the mixing matrices 
u^*(£p)u^(£ a ), u^*(£p)u^(£ a ), and u^*(£p)u^(£ a ), whose coefficients are less suppressed by a 
small x parameter, and that the rates for /1 — > eee and fie conversion in nuclei are proportional 
to \u%*(e)u%(n)\ • Since for instance in scenario A, |wf(£ a )| = \u^{l a )\Vka/y-i (attach a factor 
of c+ to y3 for scenario B), we see that the dominant terms for the LFV muon decays and fie 
conversion in nuclei are controlled by the combination £, e n = u^*(e)u^(fl) [and similarly in 
the IH case by £ ejtl = u^*(e)Ug(/i)]. We therefore seek for regions in which the combination 
£, ejl would be significantly diminished. For the mixing angles and the neutrino squared mass 
differences we use the central values from the global fit in Ref. H6l : sin 2 Q\2 = 0.32, sin 2 023 = 
0.50, Am^ = 7.6 x 10~ 5 eV 2 , {Am^ \ = 2.4 x 10~ 3 eV 2 , and set 13 at its upper limit, sin 2 6u = 
0.05. We choose OC\ — OC3 = while leaving the Dirac phase 8 and Majorana phase OC2 free- For 
the SM parameters we use the numbers of the Particle Data Group. We find that £ e[1 approaches 
its minimum at (5, atq) ~ (n/2,2n) for the NH case and at (5, Oq) ~ (0.64,3.66) for IH. This 
result is independent of the x parameter (for x not too large in scenario B). 

In Fig. [5] we show the branching ratios of the muon decays and fie conversion in nuclei and 
the largest tau decays in the NH case by scanning over the Dirac phase 8. We have assumed 
x = 8 x 10~ 3 and 0C2 = 2k. The corresponding result for the IH case is depicted in Fig. [6] as 
a function of OC2 at the same x parameter and 8 = 0.64. One sees from these two figures that 
without breaking the stringent bounds on the muon lepton the branching ratios for some of the 
leptonic tau decays can approach the level of 10 11 for almost the whole range of the scanned 
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Figure 8: Branching ratios sampled over (x, OC2) (upper panels) and (x, 5) (lower panels) for IH 
in scenario B and best-fit mixing angles. 

phase. This is much enhanced compared to the case of tribimaximal mixing shown in Figs. [3] 
and|H but is still 3 orders of magnitude below the current sensitivity. 

The above tendency encourages us to scan over a larger set of parameters. So we finally 
sample over the x parameter from 8 x 10 to 8 x 10 -3 and one of the phases in its whole range 
while keeping the other phase fixed at the value that minimizes £ eiU . Our results are shown in 
Fig. [7] for the NH and in Fig. [8] for the IH case respectively. In both figures, the upper panel 
scans over x and the Majorana phase OC2 while the lower one is over x and the Dirac phase 
8. We include only the most stringent \ie conversion on gold and the largest tau decay that is 
available in each case. In the most optimistic situation, some x decays can reach the level that 
is about 2 orders of magnitude below the current sensitivity. We notice that the two figures in 
the same panel have a similar pattern. This arises because the decay x — > ia^p^-p with a P 
is dominated by one Feynman graph which is almost the same as any of the two graphs for the 
decay x ->■ l a l a Za- 

5 Conclusion 

The origin of tiny neutrino mass has remained mysterious after years of endeavor. From the 
viewpoint of effective field theory the tiny mass can be accommodated by the canonical seesaw 
mechanisms. But it is generally hard to explore in current experiments the physics that would 
be responsible for the mechanisms because the relevant physics scale is very high and the new 
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interactions with light particles are generally too weak. It is thus highly desirable if there is any 
new mechanism that would predict accessible effects beyond the neutrino mass. 

There are two basic approaches to relax the tension between the accessibility of new physics 
and the effectiveness in producing tiny neutrino mass. One can either attribute the mass to a 
higher order quantum effect or postpone its appearance to a higher-dimensional effective in- 
teraction. An explicit model has been recently attempted in the second approach ifTOl . The 
idea is to avoid the conventional dimension-five interaction by composing new fields in higher- 
dimensional representations so that the first contribution to the neutrino mass occurs at dimen- 
sion seven. The new particles enjoy the SM gauge interactions, and thus if not very heavy 
would be produced at high energy colliders like Tevatron and LHC. The point that we want to 
emphasize here is that to establish the kinship of those particles to the origin of neutrino mass it 
would be necessary to detect their interactions with light leptons. These interactions are as usual 
shaped by the mixing effects between the light and heavy particles, and thus should also leave 
their fingerprints in precisely measured flavor-changing processes at low energy. The purpose 
of the current work has been to examine if there is any chance to look for the mixing effects in 
LFV transitions of the charged leptons. 

We have made a systematic analysis of the model. In particular, we provided a convenient 
parametrization of the leptons' mass matrices in terms of independent physical parameters. By 
diagonalizing them explicitly the lepton flavor structure becomes transparent in interactions. 
The contributions of these interactions to the radiative, purely leptonic decays of the charged 
leptons, and the fie conversion in nuclei are then computed. We considered how the stringent 
constraints from the muon lepton affect the decay processes in the tau sector. Generally the 
current experimental bounds on the decay /i — > eee and the \le conversion in nuclei, in particular 
are so strong that it is very difficult to observe the tau lepton decays. However, our sampling over 
the unknown phases and Yukawa coupling shows that there are small regions in the parameter 
space in which some tau decays have a branching ratio that is about 2 orders of magnitude 
below the current bounds. It will be challenging, if not hopeless, to observe in those decays the 
mixing effects related to the neutrino mass generation. 
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